Report 2308 


4 bo) 
SHinetoN® 


HYDROMECHANICS 
O 


AERODYNAMICS 
Oo 


STRUCTURAL 
MECHANICS 


Oo 


APPLIED 
MATHEMATICS 


oO 


ACOUSTICS AND 


ION 


2308 


x 
7 


( DOCUMENT | 


\ SN gamegiias ff 


THE PRESSURE ON THE SEA BOTTOM DUE TO A 
MOVING PRESSURE DISTRIBUTION 


by 


V.J. Monacella 


and 


J.N. Newman 


Distribution of this document is unlimited. 


HYDROMECHANICS LABORATORY 
RESEARCH AND DEVELOPMENT REPORT 


January 1967 Report 2308 


DAVID TAYLOR MODEL BASIN 
WASHINGTON, D. C. 20007 


THE PRESSURE ON THE SEA BOTTOM DUE TO A 
MOVING PRESSURE DISTRIBUTION 


by 


V.J. Monacella 
and 


J.N. Newman 


Distribution of this document is unlimited. 


January 1967 Report 2308 
S-F011 02 32 
Task 2382 


TABLE OF CONTENTS 


ABSTRACT ite icicle i) Rese NEON HT, siti scotia vi races atugeteatsiesstslcaucnsactacrnstsesate aa eae tees 


INDIO ON IESE RST UIA, VON MO) CY VANE LOIN cccacecaceqonseanccheccavaaceqaobasceo6secoossoashaabasohauasacnedenocaqzenana0364 


INTRODUGCIION isgotsudiccese eae esaceese ss tensieaie toes eclesecdeen sth at cae use atanitactasta it clesc [eh U sau cancanie a eae 


MATHEMATICAL FORMULATION OF THE PROBLEM... cere ceeeeteees 


ON EA ON CORN Oa) 0500) 12100) 8G) O71 oc ccoscasesdeseqscancno4en600005460.n00006 aq00o9500000q good HASH UdEDaEBOSAEARDOAZoARZGZAOGC 


SP SPA EDON ANTE Nee AN) 251 27 OD. INVA (OM sc cconecoaasoscoo86b0ce0sancondescnodoseadaecenenonacrnoasanodpobyoecaccqoboqN6° 


NUMBERIGAISRES UID iii ese cee cicccetec vs csccte tne essetcaatin errata ante Ssh aataaelpsesianseeetcces emeremneean 


APPENDIX A — NUMERICAL EVALUATION OF J, AND 1, oo. 


Jed) OHH OY fA) DI |( Ol Shek cor aac ohare anbssaannehcasepesanocac aocrbeackioassedensobaseobs Heaton sépsuecaqdbouaaBSsonecesoudg aduoocacodaocoa00%e> 


Figure 


Figure 


Figure 


Figure 


Figure 


Figure 


Figure 


Figure 


Th = (COG TOHENES SKVSTHEIM ‘cccrocoossxccsacdeenaboadaocooss6b6se00cs000e9s0as6s000as0700s06029cE0RaAerIoGaNCAeoONCE 


B 

2 — Bottom Pressure for — =1, F = 0.5,(7 = 0) Beene es teeseadasnacsoveadoséec 
h B 

3 — Bottom Pressure for ioe i, i? = 0.85, Fa (Oem WMP ee reenent cot eaccnoscétaddacecds 
h B 

4 — Bottom Pressure for 21s PS 1G be SS StO hl caiial, Sere cease 
h 

5 — Bottom Pressure for a. Poel iy peer Pa (=I Bra | Isa no eee eR enna neeecaosesbobo0ac0 
h B 

6 — Bottom Pressure for = 0.6, F = 0.85, (C = 03) RRR soya caodncpbene 


h B 
7 — Bottom Pressure for ee 0.6, F =1.5, fs = 0.3) BP so tonoaoaueRtat 


h B 
Speeds for rie 0.6, ae 0, (7 = 08] sneha une peste alsnenes sy saeeeeeee ame ce Ree ee 


L L 


10 


10 


11 


11 


12 


h B 
Figure 9 — Bottom Pressure for aris 0.25, F = 0.625, e = as] See aeey mae ama 
: h B 
Figure 10 — Bottom Pressure for rie O06, 17 <= Ae, aan (Caen Mie A acho eset dou Gsdndbasuaoes 


Figure 11 — Comparison of Exact Theory with Shallow-Water Approximation 


y h B 
fOr =) Oe —— OMEN OD OSA oleate AO u Mb asscsscvantuwrsaescsacm nce scee. 
L Ib L 


LIST OF TABLES 
Table 1 — Nondimensional Parameters Used in Examples .............0:ceccessesseeeseesseseseeeeees 


Table 2 — Dimensional Parameters Used in Examples ..0...0......cccccccecseeseeseessseesseesteeeneees 


ili 


NOTATION 
Beam of pressure distribution 


Speed 


Froude number ( =) 
Vgh 


Gravitational acceleration 


Fluid depth 
h 
Nondimensional depth (2 -) 


Length of pressure distribution 
Pressure 


Constant base pressure 


Rectangular coordinates 


Beam-length ratio (B/L) 


Nondimensional coordinate (2 +) 


Nondimensional coordinate (2 


Fluid mass density 


ABSTRACT 


The dynamic pressure on the sea bottom due to a constant, rectangular 
pressure distribution moving at a steady speed over the calm-water surface 
is determined. The problem is formulated and solved within the framework 
of potential-flow, linear-wave theory. Numerical results are presented for 
one beam-to-length ratio, for various water depths, and for both subcritical 


and supercritical speeds. 
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INTRODUCTION 


A pressure distribution moving over the calm-water surface will generate both a local 
disturbance around the distribution itself and a wave system which propagates behind it 
similar to those caused by a moving ship. Such a surface disturbance will obviously affect 
the pressure on the sea bottom to an extent depending on the shape and magnitude of the 
surface disturbance, its velocity, and the depth of water. 

The velocity potential associated with a pressure distribution of arbitrary shape moving 
on the surface of a fluid of finite depth is given by Wehausen and Laitone! and, by using 
Bernoulli’s equation, the pressure can be found anywhere in the fluid. If, in particular, the 
surface distribution is restricted to a rectangular area of constant pressure and if the potential 
is evaluated on the sea bottom, the analysis is somewhat simplified. Nonetheless, the pres- 
sures on the bottom are given in the form of integrals which must be numerically approximated. 
The necessary numerical analysis and computer programming have been carried out, and the 
results are presented herein. 

Insofar as a moving rectangle of constant pressure is a reasonably good approximation 
to a ground-effect machine, the bottom pressures presented herein may be interpreted as those 
due to the passage of such a vehicle. A knowledge of these pressure variations is of inter- 
est in the study of pressure-activated mines. Such interest provided the motivation for this 


investigation. 


MATHEMATICAL FORMULATION OF THE PROBLEM 


Assume that a constant pressure p = py is distributed over a rectangle of length L and 


width B and that this pressure distribution is moving over the free surface of a fluid of depth 


References are listed on page 16. 


h at a steady speed c parallel to its length. It is required to find the dynamic pressure every- 
where on the bottom and for all time. Since the motion is steady with respect to a coordinate 
system moving with the pressure distribution, the problem can be treated as being independent 
of time; thus (2, y, 2) can be taken to be a right-handed coordinate system moving in space with 
speed c in the positive # direction. Alternatively, both the pressure distribution and the co- 
ordinate system can be taken as fixed in space and with a uniform stream moving with speed 

c in the negative z direction. In either case the origin of the coordinate system is located 


in the plane of the undisturbed free surface, the 2-coordinate is taken to be positive up, and 


the pressure distribution occupies the region |a| < Re ly| < mi (Figure 1). 


For convenience the constant atmospheric pressure will be taken to be zero so that everywhere 
on the surface outside the rectangle p = 0.* 
With the usual assumptions of an inviscid, irrotational fluid there exists a velocity 


potential ®, whose positive gradient is equal to the fluid velocity vector; it can be written 


® (2, y, 2)= (a, y, 2) — cx 


where (2, y, 2) represents the perturbation potential. It will be further assumed that the 
disturbance on the surface is sufficiently small so that the problem can be linearized. 
The potential ¢ (2, y, 2) must satisfy the following conditions (see Reference 1): 
1. v* ¢ =0 everywhere in the fluid (Laplace’s equation). 
2 
@ c 
b, + — oy, = — Pp, on 2 = 0 (linearized free-surface condition). 


bo 


3. @, =0 on z =—h (rigid-bottom condition). 
Px 


4, by =0 (2 1/2) as vw > + « (radiation condition). 


?, 


where c is the forward speed of the pressure distribution (alternatively, c = speed of the 
free stream); 
for jel < — and lyl < = 
or |z = ain = 
Po = 9 Yiles 9 
pe } 
0 otherwise 


p is the fluid density; 


g is the gravitational acceleration; and the subscripts designate partial differentiation. 


*Note that there is no loss of generality provided the base pressure Po is defined relative to atmospheric 


pressure, 


(0, 0, —h) 


Figure 1 — Coordinate System 


SOLUTION OF THE PROBLEM 


The potential satisfying the conditions outlined in the previous section is given! by 


B/2 /2 1/2 
f f azay f d@ sec 0 


WIPO 8 f} fA fh, fD 0 


p (2x, Y, 2) = 


kh cosh k(2 +A) sech kh 
Cp 19 [2 (O = BD) COS 


1 2 
0 kh - — sec* @tanh kh 


F2 
cos [k (y -¥) sin 6] [1] 
B/2 L/2 7/2 
Po 
: azay f d6 sec 6 
7™pCc 
-B/2 -L/2 95 


k, cosh k, (2 + A) sech kh 
ee ee COs la, (Se) Cos O)) 


1 
1 - — sec* @sech? ko h 
F2 


cos [k, (y -7) sin @] 


c 
—— is the Froude number based on depth 


where F = 
gh 
1 
ky = kg (9 is the positive, real root of kh - — sec” @ tanh kh = 0 
F2 
-if(1 
cos —] forF>1 
6 EF 
0 
0 for F <1 


and means that the integral is to be interpreted as a Cauchy principal value. 
If z is set equal to —A, and the integration with respect to (@, 7) is carried out, then 


since 


b/e sin [k(# - 2) cos 6] 9 L 
dz = 5 sin (: Z 608 0) 


TT cos [k(x - 2) cos 6] k cos 


sin (& x cos i 


cos (k x cos @) 


and 


B/2 
2 
dy cos [k(y - y) sin 6] = sin (k B/2 sin 6) cos (k y sin 6) 
k sin 0 
-B/2 
it follows that 
4 Po 7/2 
dh (@, y,-h) = Hl d0@ sec 0 csc 0 
7™ pc 0 
” k-1 h sech kh ; 
+ dk a a sin (k L/2 cos 6) 
) kh - — sec? @ tanh kh 


F2 


sin (& B/2 sin 6) sin (& x cos @) cos (ky sin 0) 


7/2 = 
4D, f k,  sech ky kh 
- GQH'SC COCs —— —————— 
7pc iL 2 2 
—- — sec sec 
5 1 6 h kgh 
F2 


sin (k, L/2 cos 6) sin (k, B/2 sin 6) cos (k, x cos 6) cos (ky y sin 8) 


From Bernoulli’s equation, neglecting the steady hydrostatic pressure and second-order 
terms, the pressure on the bottom is given by p = pc¢,. Thus the ratio of the bottom pres- 


sure to that of the constant surface pressure p, is given in nondimensional form by 


D 4 4 
=e (Sy) = — (Ge) CS [2] 
Po m2 7 
where 
7/2 co 
H sech KH 
hell d@sec 0csc 6 dk <a ee 
0) 0 KH —- — sec? 6 tanh KH 
F2 


sin (K cos @) sin (K B sin 6) cos (K € cos 9) cos (K 7 sin 9) 


7/2 


i= [ d@sec 6 csc 0 
J, 1 - — sec* @sech? K,H 


F2 
sin (K ,cos 6) sin (K,f sin 6) sin (Ky €cos 6) cos (K47 sin 8) 


sech KH 


5 


Here B = B/L is the beam-length ratio 
& x 
Dy SLL Of 
H h 


1 
K ,H - — sec? 6 tanh K,H =0 
F2 


and kK )(@) satisfies the equation 


It should be noted that /, and /, are, respectively, even and odd functions of é and that they 
both are even functions of 7. Therefore since 


4 4 
Ee ee a) os Ga) 
Po 7 7 


it is sufficient to evaluate the integrals for only positive values of their arguments. 


SHALLOW-WATER APPROXIMATION 


The limiting value of the pressure ratio p/p, as the depth approaches zero is referred 
to here as the shallow-water approximation. From Equation [2] it can be seen that if H +0, 
KH +0, and F’ + »«, then 


7/2 


Bit’ fag ape 
Hee ag ea d@sec 0csc 0 — sin (K cos 9) sin (K B sin 6) 
Po 7 Kk 
0 0 
cos (K €cos 6) cos (K 7 sin 4) 
1 Vane = _ ~_— 
mee [p (S n) +P (-6 mn) +P (é -n)+P(-S -n)] 
where 
1/2 fos) 
cas 4 CK. ; ; 
Di(é 7) = — d@sec 6 csc 6 —— sin [K(1 - €) cos 6] sin [K(B —7) sin 6] 
2 Kk 
i 0 0 
17/2 
2 X +tan @ 
= — d@sec 6csc 61n| ————_| sgn [(1- €) (B-7)] 
2) ; X -—tan 0 


and 


1 
Here sgn (x) = +1 according as x 2), respectively. With the substitution g = x tan 0 
it follows that 


1 
4) sen ((1 - €) (8 - 7) = sen [1 - €) (8 - 0) 


This leads to the simple result that 


Bas hoon || <al and In| <B 
1 if |é|>1 and/or |n| > B 


that is, as one would expect, the pressure p is zero everywhere except directly under the rec- 


tangle where it takes the value p,. 


NUMERICAL RESULTS 


As stated before the integrals presented in Equation [2] can only be evaluated 
numerically. A program for this purpose was written for the IBM 7090 and several examples 
are given graphically in Figures 2 through 11. Details of the numerical analysis are outlined 
in the Appendix. All of the calculations are for a beam-length ratio (8) of 0.3. Table 1 
shows the values of the remaining parameters pertinent to the various figures. 

In order to give some dimensional concept to the figures one can assume the rectangle 
to be 100 feet long and 30 feet wide. With this configuration, Table 2 shows the dimensions 
equivalent to the parameters of Table 1. 

For subcritical speed* (F < 1.0) there is both a transverse and divergent wave system 
propagated behind the rectangle. As would be expected the oscillations of the centerline 
(y/L = 0) pressure traces of Figures 2, 3, 6, and 9 are similar to those of the transverse 
wave system. For supercritical speeds there are no transverse waves and, consequently, 
one would not anticipate an oscillatory behavior of the centerline trace far downstream. 
Figure 8 shows a comparison of centerline pressures for subcritical and supercritical speeds. 
Off the centerline, however, the pressure trace can still exhibit an oscillatory pattern due 
to the diverging wave system. Since the diverging waves are, in general, of shorter length 


than the transverse waves, the extent to which they affect the bottom pressure depends very 


c 


Tee 


*Critical speed = 


largely on the depth of water. As can be seen from the figures, typical behavior of a super- 
critical speed, centerline pressure trace is a building up (from zero) of the pressure near the 
forepart of the rectangle, a falling off to a negative peak behind it, becoming positive again, 
and finally approaching zero asymptotically. 

The data presented in Figure 11 is intended as a comparison with the shallow-water 
approximation. In this example y/L = 0, while the speed is held constant. The traces are 
for Froude numbers of 3, 5, 7, and 10. (Recall that in the shallow-water approximation H > 0 
while F + ~.) As F increases the abrupt pressure jump at the bow and stern becomes more 
evident, as predicted in the shallow-water approximation. Also with increasing F the bottom 


pressure directly under the rectangle appears to be approaching the limiting value of one. 


TABLE 1 TABLE 2 


Nondimensional Parameters Used in Examples Dimensional Parameters Used in Examples 


Distance Abeam 
of Centerline 


ft 


Figure Depth Speed 


ft knots 


0, 1.0, 2.0 


9 0, 100 
0, 1.0, 2.0 ; 0, 100, 200 
0, 1.0, 2.0 ; 0, 100, 200 
0, 1.0, 2.0 5 0, 100, 200 
0, 1.0, 2.0 ; 60 0, 100, 200 
0, 0.5, 1.0 ; 60 0, 50, 100 
0, 0.25, 0.625 ; 5 0, 25, 62.5 
0, 0.25, 0.5 " 05 0, 25950 
11 | 10.0, 20.4, 0 
40.0, 111.1 


NOTE: L = 100 feet and B = 30 feet. 


DIMENSIONAL PARAMETERS 
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2 Ht 
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~0.04 yt 1 r 
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x | 
~0.06 laa 
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a h B 
Figure 2 _ Bottom Pressure for a 1.0, F = 0.625 (2 = 0.3) 
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Figure 3 — Bottom Pressure for ian IFO FE = 0585 = 0.3) 
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DEPTH = 100 FEET 
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Figure 4 — Bottom Pressure for — =1.0, F =1.5 (- = 
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PARAMETERS 


BEAN 30 FEET 
SPEED = 50.4 KNOTS 


DISTANCE ABEAM 
OF CENTERLINE IN FEET 


0 4.0 5.0 


0.3 
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Figure 5 — Bottom Pressure for ile 1.0, F =2.0 z =0.3 
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Figure 6 — Bottom Pressure for ait 0.6, F = 0.85 rs = 0.3 


PRESSURE 
DISTRIBUTION 


1 ' 
ee oe DIMENSIONAL PARAMETERS 
LENGTH = 100 FEET BEAM = 30 FEET 


DEPTH = 60 FEET SPEED = 39.0 KNOTS 


DISTANCE ABEAM 
OF CENTERLINE IN FEET 


0 
a5) 
——-—100 
4 i a 
1.0 2.0 3.0 4.0 5.0 


11 


DISTRIBUTION 


-5.0 


WS? (Ga ae Sali Map al =] 


-4.0 -3.0 -2.0 -1.0 0.0 1.0 


DIMENSIONAL PARAMETERS 


CENTERLINE PRESSURF 
LENGTH = 100 FEET BEAM = 30 FEET 
DEPTH = 60 FEET 


SPEED IN 


Figure 8 — Comparison of Bottom Pressures at Subcritical and Supercritical Speeds 
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Figure 11 — Comparison of Exact Theory with Shallow-Water Approximation 


APPENDIX A 


NUMERICAL EVALUATION OF 7, AND /, 


With the exception of special cases, such as the shallow water approximation, the 
integrals /, and /, must be evaluated numerically. A FORTRAN program was written for this 
purpose, and a summary of the numerical analysis is presented here. For simplicity, and 
where the meaning is clear, the integrands of the integrals that follow will be omitted. 


It was pointed out before that K (6) is the positive, real root of the transcendental 


1 : 
equation KH ~ — sec” 6 tanh KH =0. This root is obtained for each required value of 6 
F2 
by means of the Newton-Raphson method. ? 
Although the analysis pertaining to the integrals is similar in both cases it is some- 


what more convenient to treat subcritical and supercritical speeds separately. 


SUBCRITICAL SPEED (F <1) 


This is the simpler of the two cases, primarily because 6) = 0. For this reason together 
with the fact that the integrands of both /, and /, are functions of cos 6, the @ integrations 


lend themselves to the Tchebysheff quadrature? in the form, 


7/2 n 


40 f(cos )= — Df (cos 4) 
n 


0 j=l 


27-1 
4n 
ture, it is especially advantageous for computer usage in that both the abscissae and (con- 


where 0; = 


7. In addition to the good accuracy associated with the Gaussian quadra- 


stant) weight functions can be readily derived within the program. It is used initially for /, 
and the @ integral of /, by setting n equal to 10. Then the integrals are reapproximated with 
n doubled. If the estimates are not sufficiently close as determined by a convergence 
criterion, n is redoubled, and the process is continued until convergence is obtained. 

For each value Fe of the double integral /, the root kK, is found and the integration on 


K is performed by first writing the infinite integral as the sum 


2K, 2K y+ 2 2Ky+4 
f dk + dk + dK + «+s 
(0) 2Ky 2K +2 


and approximating each of the integrals in the series by a repeated application of the 10-point 


Gauss-Legendre quadrature;* that is, the range of integration of each integral is repeatedly 


14 


(subdivided, and the quadrature is applied to each subdivision until convergence is obtained. 
The series (of integrals) is terminated when the ratio of the absolute value of any one integral 
to the absolute value of the sum of terms up to that integral is sufficiently small. The first 
term in the series involves the calculation of a principal value integral but since an even 
order (10 point) quadrature is used in which the abscissae are symmetric about (but do not 
include) K, the singularity there can be ignored. * 


SUPERCRITICAL SPEED (F > 1) 


1 
The significant differences when F > 1 are that 0) = cos” ! (=) , and there is a log- 
arithmic singularity at @ = 6, in the integral /,. Because of these the Tchebysheff quadrature 
is not applicable and the 10-point Gauss-Legendre quadrature is used. 
The integral /, is approximated by the same method described previously, namely, sub- 
dividing the range (9), 7/2) of integration and applying the quadrature over each sub- 


division, the number of such subdivisions being whatever is necessary to obtain convergence. 
The integral /, is written as the sum 


lh Siig tling 


where 
95 oo 
lo = j dé | dk, and 
0) ty) 
7/2 00 
Ina i dé f dK 
% ) 


In both /,, and /,, the 6 integral is approximated by using the Gauss-Legendre quadrature in 
combination with the 10-point Gaussian quadrature for integrals with a logarithmic singularity. > 
The integral is approximated in the same way as /,, except that the Gauss-Legendre quadra- 
ture is used over all the subintervals except the one including 6, for which the logarithm 
quadrature is used. 

The K integral of /,, is calculated exactly as in the subcritical case. However, for 
I,, it can be readily seen that there is no non-zero root of the transcendental equation and the 
integrand is well behaved. For this case the integral is written in the same infinite series 
as before and is approximated exactly the way as for subcritical speed with the exception 


that Ky =0. 


By the requirement that sufficiently close estimates or of convergence be obtained, a 
relative error of less than 1 percent is typically implied. 
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